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The phenomena of radial and axial segregation in a horizontal rotating cylinder containing a 
mixture of granular particles of two different species have been modeled using discrete particle 
simulation. Space-time plots and detailed imagery provide a comprehensive description of what 
occurs in the system. As is the case experimentally, the nature of the segregation depends on 
the parameters defining the problem; the radial component of the segregation may be transient or 
long-lasting, and the axial component may or may not develop. Simulations displaying the different 
kinds of behavior are described and the particle dynamics associated with the axially segregated 
state examined. The importance of an appropriate choice of interaction for representing the effective 
friction force is demonstrated. 

PACS numbers: 45.70.Mg, 45.70.Qj, 64.75.-l-g, 02.70.Ns 
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I. INTRODUCTION 

The mechanisms responsible for some of the more sur- 
prising features of granular matter continue to challenge 
understanding. In view of the fact that the behavior en- 
countered often differs strongly from systems governed 
by statistical mechanics and thermodynamics, there is 
little by way of intuitive help from such theories in trying 
to understand granular matter, both at rest and in mo- 
tion. Segregation is perhaps the most prominent of these 
features; the ability of noncohesive granular mixtures to 
segregate into individual species, despite the absence of 
any obvious cntropic or energetic benefit, makes this a 
particularly fascinating phenomenon. 

Granular separation and mixing are central to many 
kinds of industrial processing that span a broad range of 
size scales; the capability for either causing or preventing 
segregation can be central to the viability of a particular 
process. Substantial economic benefits could be gained 
from a systematic understanding of the complexities of 
granular flow at a level similar to fluid dynamics; the 
outcome would be improved control, reliability and effi- 
ciency, a notable advance over currently used approaches 
that are often empirical. Analogous segregation processes 
occur in nature, where an understanding of the under- 
lying mechanisms might benefit resource management. 
Consequently, considerable effort has been invested in 
exploring these phenomena, but, absent a general the- 
ory of granular matter, much of the progress in this field 
[Hi S [Sj 4] relies on computer simulation. 

There are a variety of conditions under which granular 
mixtures segregate; these include shear @ , vibration 
0, H, and rotation [13, [O, [ll- Experiments on ro- 
tational segregation, one of the more extensively studied 
systems experimentally and the subject of the present pa- 
per, employ a horizontal cylinder revolving at constant 
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rate and partially filled with a mixture of two species of 
granular particles. Given the appropriate kinds of parti- 
cles, and under suitable conditions, the mixture sponta- 
neously segregates into axial bands of alternating com- 
position. A substantial body of experimental data has 
been assembled, but because of the numerous parame- 
ters entering the experiment only a limited portion of 
the multidimensional phase diagram has been explored. 
The relevant parameters include the nature of the gran- 
ular particles themselves, ranging from the most familiar 
forms of granular matter such as sand, through various 
organic products, plastic and glass beads, to metal ball 
bearings. The particles have a broad spectrum of elastic 
and frictional characteristics, and sizes that range over 
several orders of magnitude. The material can be in a 
dry or wet state, the latter case known as a slurry. The 
cylinder can have different shapes, most commonly one 
with a uniform circular cross section. Further parameters 
include the cylinder diameter and length, the rotation 
rate, the relative fractions and sizes of the granular com- 
ponents, and the filling level of the container. In marked 
distinction to fluid systems, there are no known scaling 
relations that allow this parameter space to be reduced, 
and while some systematic trends have been noted, there 
is little guidance available for predicting the outcome of 
an experiment involving some unexplored combination of 
parameters. 

The most direct observations address the structure of 
the upper surface while more elaborate studies ex- 
amine the interior using, for ex amp le, noninvasive mag- 
netic resonance imaging (MRI) [ij, [3, [l^ . Axial segre- 
gation may be complete, with each band composed ex- 
clusively of one particle species, or it may be only par- 
tial, with bands characterized by higher concentrations 
of alternating species. There is an additional effect that 
occurs entirely beneath the surface leaving little exter- 
nal evidence, namely radial segregation, which produces 
a core rich in small particles extending along the length 
of the cylinder. Radial segregation occurs early in exper- 
iments starting from a uniformly mixed state, and may 
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or may not persist once axial segregation appears. Axial 
band formation might be regarded as a consequence of 
this core intermittently bulging by an amount sufficient 
to reach the outer surface |f4| . The axial band patterns 
can also exhibit time dependence; examples of such be- 
havior involve coarsening, in which narrower bands grad- 
ually merge into broader bands [16], and the appearance 
of traveling wave patterns on the surface JT,, J8| . Fur- 
ther discussion of the experimental history appears in 
[13, 0) [3 ' tiut information about this phenomenon still 
continues to accumulate. More recent results showing 
the richness of the segregation effect and its dynamics in 
dry and/or wet mixtures (slurries exhibit behavior that 
is similar in many respects, despite the lubricating role 
of the fluid medium) are described in [IS, HHj Ull • The 
problem can be extended in different directions, for ex- 
ample (23j considers mixtures with three and more gran- 
ular components. That there is little opportunity for 
scaling to play a role in understanding this phenomenon 
is demonstrated in [23|, where the ratio of the cylinder 
and particle diameters determines whether axial segrega- 
tion occurs and if its appearance depends reversibly on 
rotation rate; the fact that the wavelength of the axial 
pattern does not scale with cylinder diameter is described 
in m. 

Given the complexity of the experimental situation and 
the many lacunae in parameter space, it is hardly surpris- 
ing that theoretical progress has been limited. Continu- 
ous one-dimensional models have been developed [26] in 
which the dynamical variables are the local concentration 
difference and the slope of the free surface; these describe 
the early stages of segregation with traveling bands and 
subsequent band merging, although experiment (27j has 
questioned their ability to characterize the process cor- 
rectly, and the obviously important three-dimensional as- 
pects of the phenomenon are absent. The problem has 
also been studied using a cellular automata approach 
based on a highly simplified model [1^. 

The most detailed theoretical approach involves direct 
modeling of the dynamics of granular particles, employ- 
ing the same computational methods used in molecular 
dynamics simulation of atomistic systems. While there 
have been a considerable number of such studies, cov- 
ering a variety of granular phenomena, there have been 
very few simulations of granular particles in a revolving 
cylinder exhibiting axial segregation. The first of these 
was an early, rather limited treatment |29| . This was fol- 
lowed by a broader study [l^], that is extended by the 
present work, in which the formation of multiple axial 
bands was demonstrated, and various aspects of the seg- 
regation phenomenon, such as band merging, the absence 
of mixing in presegregated systems, and the manner in 
which the behavior depended on the choice of parame- 
ters, were investigated. The model differed from that of 
the present paper in that the transverse restoring force 
(discussed in detail below) was absent and a Lennard- 
Jones type of normal repulsion was used rather than a 
linear force. Exactly the same kind of axial band for- 



mation was subsequently reported in [30| (although this 
similarity was not explicitly noted) using the same kind 
of model, but with the linear force. The appearance of 
radial segregation is described in ^Sl*] , but only for a two- 
dimensional system. 

The purpose of the present paper is to analyze a 
discrete-particle model that exhibits both axial and ra- 
dial segregation, in which the appearance of the radial 
effect precedes the axial, and then either disappears or 
persists as an axial core. Experimentally, it is the smaller 
particles that are found to congregate in the interior, and 
this feature is reproduced by the simulatioiis; one issue 
that was not resolved in the earlier study [I4] was the 
fact that the opposite behavior was encountered, but it 
will be shown that this is due to the choice of interparti- 
cle forces. The goal here is not a systematic coverage of 
parameter space, since the computational resources for 
this would be considerable, but rather a demonstration 
of the typical kinds of outcome with just a few examples 
of parameter dependence; indeed, due to a lack of de- 
tailed reproducibility under different initial conditions, as 
will become apparent subsequently, obtaining a compre- 
hensive picture of the "average" behavior would require 
multiple simulation runs for each parameter set. Typi- 
cally, the majority of runs described here involve larger 
systems than before, with a bigger particle size differ- 
ence and a more slowly rotating cylinder; these changes, 
together with the additional work required by the force 
computations, called for a greater computational effort 
than previously. 



II. SIMULATION METHODOLOGY 
A. Granular models 

The models used for granular simulation are generally 
based on spherical particles whose resistance to overlap 
is expressed in terms of a continuous potential function 
HjIH,!!!]. In addition to this excluded- volume repulsion, 
particles are subject to forces that aim to replicate the 
effects of inelasticity and friction. There are limitations 
to the accuracy of such models, with one empirical mea- 
sure of success being the degree to which the essence of 
the phenomenon under study is captured. 

Consider a pair of granular particles i and j with diam- 
eters di and dj , respectively. The repulsive force between 
particles depends linearly on their overlap, 

fv — {^ij ) '^ij : ^ij ^ d^j ^ (1) 

where = — Tj is the particle separation, = 
and dij — {di + dj)/2 is the effective diameter. An alter- 
native is the Hertz interaction that depends on the 3/2- 
power of the overlap [11], whereas in [l^] the repulsive 
part of the Lennard-Jones interaction was used. Dissi- 
pative forces act for the duration of each collision. The 
first of these is velocity-dependent damping along the line 
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between particle centers, 



(2) 



that depends on the component of relative velocity Vij = 
Vi — Vj parallel to . Here, 7„ is the normal damping 
coefficient, assumed to be the same for all particles. The 
total normal force between particles is then /„ = fv + fd- 
Frictional damping acts transversely at the point of 
contact to oppose sliding while particles are within inter- 
action range. 



/. = -min(7r^|<.U=-^^|/„|)*,^ 



(3) 



where the relative transverse velocity of the particle sur- 
faces at this point, allowing for particle rotation, is 



(n. 



(4) 



and (jJi is the particle angular velocity. The value of 7s' ^ 
depends on the particle types and c_,-, and is the 
static friction coefficient that limits the transverse force 
to a value dependent on | /« | . 

In a model of this type there is no true static friction, 
a practical consequence of which would be stick-slip mo- 
tion. A way of at least partially overcoming this limita- 
tion, though not strictly a correct means of incorporating 
the effect, is to introduce a tangential restoring force that 
acts during the collision and depends on the cumulative 
relative displacement at the point of contact 
This force has the form fg = —kgUij, where 



(coll) 



,(T)dT 
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is evaluated as a sum of vector displacements over the 
interval the particles have been in contact; the magni- 
tude of fg is also limited by /^^'^^|/n|- In addition, to 
avoid occasional unrealistically large displacements (al- 
though this does not appear to affect the behavior), if 
I My I > 0.1, an arbitrarily chosen limit, it is reset to zero 
(a form of ratcheting intuitively motivated by the asper- 
ities responsible for friction); an alternative treatment 
appears in ISSl . This history-dependent force was not 
present in [Tj]; the results below suggest that it plays 
an important role in achieving the correct form of radial 
segregation. The total transverse force is ft = fs + fg- 

The curved cylinder wall and the flat end caps are 
treated as rough and smooth boundaries respectively. 
The interaction of particles with the curved boundary 
(together with gravity) drives the system, so this force 
includes the same components as the interparticle force. 
On the other hand, to minimize spurious effects associ- 
ated with the end caps, only fy and fd act there. Further 
details concer ning the interactions (with the exception of 
fg) appear in [12|, together with a discussion of the fric- 
tion coefficients and the particle- wall force computations. 
Other aspects of the simulation follow standard molecu- 
lar dynamics procedures [36.] : neighbor lists are used to 



efficiently organize the force computations, the transla- 
tional and rotational equations of motion are integrated 
with the leapfrog method, and parallel processing can 
improve performance. 



B. Parameters 

Those parameters that also appeared in [l^ have 
been assigned similar values here. The first of these 
is the gravitational acceleration, g = 5, which then 
relates the dimensionless MD units used in the sim- 
ulation to the corresponding physical units. Thus, if 
Lmd is the length unit (in mm), then the time unit 
is Tmd ~ 10^-^\/5Lmi) s. If the cylinder rotates with 
angular velocity il (MD units) the actual rotation rate 
is rt/{2TrTMD) ~ 7.1fl/^/LMD Hz; for 3 mm particles, 

— 0.1 is equivalent to 25rpm, a typical experimental 
value. 

Axial segregation tends to develop slowly, leading to 
long simulation runs. Computational cost can be re- 
duced by using the largest possible integration time step; 
the step size 6t is limited by the highest particle speeds 
encountered and here the value used is St = 5 x 10~^ 
(MD units). Numerical stability considerations then set 
a lower bound on the collision time to ensure it remains 
much larger than 5t (typically 30-40x). Collisions be- 
tween particles in the slowly varying bulk interior can be 
protracted events, but the opposite is true for particles 
bouncing rapidly along the upper free surface, and this in 
turn constrains coefficients such as particle stiffness. For 
the present work, most runs use fc„ = 1000 in Eq. ([1]). 
Increasing this to, e.g., 10^ would require a 10- fold re- 
duction in St and a much longer computation; however, 
since there are ~ 10^ steps/revolution, some runs already 
exceed lO*" steps, so this is not presently feasible in gen- 
eral. 

The relatively small cylinder diameters and fill levels 
used here (details below) limit the compression forces 
that particles can experience. At the base of a static 
column of 10 small particles the overlap (which can be 
regarded as compression) will be 5% (for fc„ and g as 
given), an amount that is unlikely to alter the qualita- 
tive behavior; indeed, tangential forces will be enhanced 
for softer particles that experience a greater number of 



longer-lasting contacts. (Particles with the 



repul- 



sion that is part of the force used in [12] have an overlap 
three times larger; the ultimately greater resistance of 
this force to compression asserts itself only at five times 
the depth.) Systems with substantially thicker layers will 
require a stiffer volume interaction (and consequently a 
smaller St). 

The values of the remaining parameters from [T^ are 
as follows. In Eq. 7„ — 5. The coefficients in Eq. ^ 
are 7^'' = 10, and in general 7I'' — j^'^ = 2, except for 
a few cases where all 7^ = 10; the particle- wall values 
are identical. The relative values of the static friction 
coefficients, e.g., fj,''^ / fi'^'^ , are set equal to the ratio of 
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the corresponding 7^ values, with the larger of the pair 
equal to 0.5. Finally, the value of the new parameter 
introduced with Eq. ^ is kg — 500, although, as with 
most parameter settings, there is considerable latitude, 
and all that is required at this exploratory stage is for 
the results to appear qualitatively reasonable. 

The nominal particle diameter is the interaction cutoff 
in Eg. (01). For small particles, to maintain consistency 
with [ij 4 = 2^/6 « 1.122 (MD units), while for big 
particles db ~ bdg. The actual particle sizes are uni- 
formly distributed over a narrow range [ds — 0.2, ds] for 
small particles, and likewise for big; the mean diameter 
of the small particles is then close to unity. Particles all 
have the same density, so the big to small mass ratio (be- 
fore allowing for the random size distribution) is . The 
relative population of big and small particles is chosen to 
give equal volume fractions; for most runs b = 1.8, with 
a big particle fraction of 0.15. 

The cylinder diameter D and length L have value 
ranges 30-40 and 120-360 (MD units) respectively. (A 
consequence of the use of soft potentials is that the effec- 
tive D and L values are reduced by approximately one 
unit.) The aspect ratio L/D, a quantity having some 
influence on the number of axial segregation bands [2^, 
lies between 4 and 12. The cylinders used here are rela- 
tively narrow, with D less than 30 times the mean par- 
ticle diameter, a limitation that has been found experi- 
mentally to influence behavior [23|; even so, the number 
of particles in a simulation can exceed 70 000. Another 
connection between experiment and simulation is the di- 
mensionless Froude number Fr = Cl"D/2g, the ratio of 
centrifugal to gravitational acceleration. Experimentally 
Fr « 10~'^-10~^, and to ensure the simulations remain 
in the correct regime, the value should not be allowed to 
become too large; this sets an upper limit to f2, and for 
the simulations Fr < 0.2. 



C. Measuring segregation 

Quantitative measures of the overall intensity of axial 
and radial segregation are used to augment direct obser- 
vation. The time-dependent axial segregation, Sa{t), is 
evaluated from binned counts along the axis weighted by 
particle mass, 



used. 
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Sa 



{Ib^nbjq) -nsiq)\)q 
b^Ub -\- Us 
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where nb{q) and ns{q) are the number of big and small 
particles in a slice (width w 2.5) centered at q along the 
axis, (. . .)q denotes an average over all slices, and rib and 
Us are the totals. Binned counts are less suitable for eval- 
uating the radial segregation, Sr(t), because, if done in 
the same way as the space-time plots discussed below, 
not all particles would be able to contribute. An alterna- 
tive measure, based on the mean-square radial distance of 
each type of particle from its center of mass, is therefore 



where (rf,) and (rs) denote the mean distances of big and 
small particles from the axis, and (r) the overall mean. 
Time is expressed in terms of the number of cylinder 
revolutions, n/j. 



III. RESULTS 

Due to the very nature of the segregation effect, not 
all of its features are readily quantifiable. The functions 
Sa{t) and Sr{t)^ Eqs. ^E^, provide a global summary of 
pattern development, but lack detailed information con- 
cerning the number of axial bands, the sharpness of the 
boundaries between segregated regions and the regular- 
ity of the band pattern. The use of space-time plots, 
coded using either color or grayscale level, provides a 
much more detailed picture of how both axial and ra- 
dial segregation evolve with time. Such plots show the 
relative populations, weighted by particle volumes, in ap- 
propriately oriented slices. For the axial plots, slices are 
normal to the cylinder axis; while this does not provide 
the same information as the surface populations seen in 
experiment, it will be similar provided that radial inho- 
mogeneity is weak. For the radial plots, slices are parallel 
to the upper free surface of the material; the nominal sur- 
face slope is determined by a linear fit to the inner 2 /3 of 
the surface away from the curved boundary, and only the 
particle populations in a slab of similar width normal to 
the surface are counted, to avoid bias due to the shape 
of the region. 

Even more detailed information can be extracted from 
snapshots of the entire system, in which particular sub- 
sets of particles can be selected for viewing; images of 
this kind can, for example, reveal interior organization 
analogous to that observed experimentally with MRI. An 
even richer visual approach employs animated recordings 
showing the full temporal development of the system. An 
animation of this type consists of an extended series of 
snapshots, recorded at regular intervals throughout the 
run; all the information needed for the analysis that fol- 
lows is in fact obtained from configurations reconstructed 
from such recordings. 

The runs described here are cataloged in Table I; they 
are labeled alphabetically for reference and listed in the 
order they first appear in the discussion. Run lengths 
were not specified in advance; a run was generally al- 
lowed to continue for as long as something interesting 
was happening, or until it seemed that the system had 
stabilized (the possibility of premature termination can 
never be ruled out), or until changes appeared to be oc- 
curring too slowly to warrant continuation. 

Most runs are begun from a uniformly mixed initial 
state, and in a few cases, from a state that is preseg- 
regated into two axial bands. The particles are placed 
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TABLE I: Summary of runs discussed in the text; the top row 
includes the defauh settings used unless indicated otherwise. 



Id° Size, otc' Forces'* R" Outcome 





L 


D 


p 


N 


b 


n 






kg 






riR^ 


Rd« 


Ax*" 


A 


240 


40 


.30 


71 600 


1.8 


.2 


M 


1000 


500 


2 




1620 


P 


12 


B 


360 


30 




57 600 




.1 












2630 


T 


17^1 


C 


160 


40 




47 300 




.1 












1800 


P 


7 


D 


120 


30 


.50 


35 900 
















470 


P 




E 


120 


30 




18 800 
















2940 


T 


6 


F 


120 


30 




18 800 














R 


1310 


T 


5 


G 


120 


30 




18 800 














R 


1310 


T 


7^4 


H 


120 


30 


.40 


24 200 
















730 


P 




I 


120 


30 


.45 


29 300 
















590 


P 




J 


120 


30 




18 800 




.1 












4290 


T 


6->3 


K 


360 


30 




57 600 
















2960 


T 


17^1 


L 


160 


40 




47 300 
















3320 


P 


7 


M 


120 


30 




18 800 










Q 






2040 


J 


3 


N 


120 


30 




18 800 















R 


1030 


I 


5 


O 


120 


30 


.45 


29 300 

















830 


P 




P 
Q 


120 
120 


30 
30 




19 200 
19 200 


1.5 
1.3 














6350 
2760 




7^5 
6 


R, 


120 


30 




18 800 








2000 








2610 


T 


6->3 


S 


120 


30 




18 800 












10 




6350 


P 


2' 


T 


180 


30 




16 800 






S 










2840 




2 


U 


180 


30 




16 800 






S 






10 




2440 


P 




V 


160 


40 


.18 


27 900 






S 










6110 




2-^4 



"Runs are denoted in the text as #A, etc. 

'Cylinder length L and diameter D, filling density p, number of 
particles N (rounded to nearest 100), big particle size 6, angular 
velocity Q. 

■^M/S: Mixed or segregated initial state. 

''Values of force coefficients: kn, kg and 7!"; all others are con- 
stant with values listed in the text. 
'^R denotes repeat of preceding run with different initial state. 
■fRun length, in cylinder revolutions. 

fRadial segregation: permanent, transient, inverted, or none. 
''Axial segregation: initial and final band count, or none. 
'Partial segregation only. 

on a lattice and assigned small random velocities (details 
that vanish after just a few collisions); particle species is 
either randomly assigned according to relative concentra- 
tion, or, in the presegregated case, determined by axial 
position. A change of random number seed allows re- 
peated runs of systems that are otherwise identical. 

In all cases the cylinder is capped at the ends. Use of 
an axially periodic cylinder leads to similar pattern devel- 
opment, establishing that cap s are not required for seg- 
regation. However, as in |12l |. the entire band structure 
is then subject to axial drift. Since there is no evidence 
that the caps affect the results (apart from relaxing the 
requirement for an even number of bands in the periodic 
case) they can be used without concern. 

In view of the extensive simulations involved, it is inter- 
esting to consider the computational performance. For a 
system of 47000 particles, and a partly parallelized com- 
putation running on a dual 3.6 GHz Intel Xeon worksta- 
tion, the simulation proceeds at a rate of approximately 
23 000 time steps/hr, equivalent to 3.7 revolutions/hr at 
n = 0.2; some of the runs therefore extend over several 




FIG. 1: (Color online) Axial and radial space-time plots for 
run #A (ur = 1620, L = 240, L/D = 6); red and blue (or 
medium and dark gray) denote higher volume fractions of big 
or small particles. Time, expressed as cylinder revolutions, is 
along the horizontal axis, and the vertical axis measures axial 
or radial (the latter normal to the free surface) position. 

weeks. Contrast this with a typical experiment that can 
be run at, e.g., 1000 revolutions/hr, irrespective of the 
number of particles. 

A. Axial and radial segregation 

The first two runs discussed are among the more ex- 
tensive carried out. While similar in respect to the ini- 
tial radial segregation that peaks after some 10-20 rev- 
olutions and the appearance of multiple axial bands, in 
one instance radial segregation persists while in the other 
it practically vanishes. The details of these two runs, 
and #5, will be described using axial and radial 
space-time plots, together with pictures of the final states 
that reveal the nature of the internal structure, and with 
graphs of segregation as measured by Sa and Sr- Both 
runs suggest that a steady final state has been achieved, 
although, given the experimentally observed long-term 
slowdown in the bond merging rate, there is no way of 
completely excluding future changes of this type, no mat- 
ter how long the run. 

Fig. [1] shows the axial and radial space-time plots for 
run ^A. Radial segregation appears very early in the 
run and a core of small particles persists throughout, 
although the outer layer rich in big particles disappears 
from the radial plot once axial banding begins; after the 
12 axial bands have formed the pattern appears stable, 
with no hint of any future change. 

Images of the final state of run ^A appear in Fig. [2l 
The first is an oblique view of the full system showing 
the 12 alternating bands. The others are views looking 
down in a direction normal to the surface that show big 
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FIG. 2: (Color online) Final state of run H^A; the full system, 
and views showing just the big and small particles (colored 
silver and gold). 



and small particles separately. These interior views re- 
veal that while the bands of big particles are not entirely 
separate, there are visible divisions between them, and 
that the bands of small particles are joined by a central 
core extending the length of the system. 

Fig. [3] shows the space-time plots for run jj^B. Here, 
unlike ^A, radial segregation is a transient effect last- 
ing less than 200 revolutions; furthermore, the early ax- 
ial pattern is not maintained, and the initial 17 bands 
are eventually reduced to 11 due to the vanishing of the 
small particle bands at the ends and the merging of two 
pairs of big particle bands (eliminating two more small 
particle bands). During the latter portion of the run, ex- 
tending over more than half its total length, there is no 
suggestion of further pattern change. Band coarsening is 
a well-known experimental result, with the band count 
falling roughly logarithmically with time [lO], although 
the dependence can be more complicated than this 22|. 

The images in Fig. [3] show the final state of ^B. Band 
separation is practically complete in this case, with es- 
sentially no misplaced big particles, and only the faintest 
remnant of a small particle core. 

The development of axial and radial segregation, Sa 
and Sr, for these two runs is shown in Fig.O The graphs 
reflect what has already been noted in the space-time 
plots and images, namely, that a preference for small 
particles in the core persists over the duration of but 
not #_B; the presence of this core affects the magnitude 
of Sa, which is smaller in #^ than in #5. 

The images in the next two flgurcs provide examples 
of radial segregation in cases where the effect is well de- 
veloped. Fig. [5] shows #C after 60 revolutions; this run 
eventually develops axial segregation. The complete sys- 
tem is shown, together with three narrow slices (of thick- 
ness 0.03L) at the midpoint and at a distance O.li from 
either end. 

Fig. [7] shows #D after 260 revolutions (graphs of Sa 
and Sr for this run appear later). The view is along 




FIG. 3: (Color online) Axial and radial space-time plots for 
#B {uR = 2630, L = 360, L/D = 12). 
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FIG. 4: (Color online) Final state of #B. 



the axis, after slices of thickness O.IL are removed from 
each end; on the left all particles are shown, while on 
the right just the big particles appear. There is a visible 
opening extending along the entire length (0.8i) near the 
center of the small particle core; while the simulations 
do not produce a totally pure core, radial segregation is 
strong despite misplaced big particles; this rather short 
run showed no hint of imminent axial segregation. 

Two general observations concerning segregation 
emerge from these (and subsequent) simulation runs. 
When radial segregation occurs it takes the form of a 
core of small particles surrounded by big particles (there 
is an exception that is discussed below); this is also the 
case experimentally, although in the simulations the core 
boundaries are not as sharp and the effect appears to 
develop more slowly. There is no preferred particle type 
near the end caps nor a tendency for bands to nucle- 
ate there; bands can form at different axial positions at 
different times; sometimes bond formation is almost si- 
multaneous and sometimes not. 
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FIG. 6: (Color online) Views of #C after 60 revolutions; the 
full system and three slices are shown. 



B. Reproducibility 

The issue of reproducibility is of particular importance 
since conclusions must take into account the variability 
of the behavior, such as different intermediate or termi- 
nal band counts. Multiple runs would be required, as is 
normal experimentally. In most cases this condition has 




FIG. 7: (Color online) Views of after 260 revolutions (see 
text for explanation). 



not been fulfilled, but an example of three runs involving 
systems that are identical, apart from the initial random 
state, is considered here (and one further example ap- 
pears later). 

Fig. [5] shows the axial and radial space-time plots for 
run together with the axial plots for the shorter runs 
#F and #G. The behavior differs (although Sa and Sr 
are only weakly affected); each exhibits early radial seg- 
regation, but they differ in regard to axial band develop- 
ment. This is just what happens experimentally under 
appropriate conditions |20|] : some aspects of the segre- 
gation {Sa and Sr) are reasonably robust, while others 
(the band details) may vary significantly between runs. 
A preferred wavelength governing the axial band pattern 
could exist, but since development history can influence 
outcome, it is probably only meaningful to refer to a 
mean number of bands (and a mean wavelength). 

C. Parameter dependence 

When comparing simulation to experiment, the param- 
eters defining the system fall into two categories. There 
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FIG. 9; Axial and radial segregation (the graphs for the latter 
peak early in each case) for runs ^E, #H, and with 
increasing fill level. 



are those that are readily changed in an experiment, such 
as the cylinder rotation rate, the fill level, and even the 
cylinder size; experimental results may cover a range of 
these parameters in order to demonstrate trends or con- 
struct phase diagrams. The others are related to the 
properties of the particles themselves; while there are 
numerous types of granular material available for exper- 
imental study, efforts have tended to focus on a rather 
limited subset, so that a systematic coverage of behavior 
in terms of material properties is presently unavailable. 

Changing the fill level can alter the nature of the seg- 
regation, as seen experimentally with slurries [2l|. Fig.[S] 
shows Sa and Sr for four runs differing only in their fill 
level, ^H, #/ and #D (because of the slight reduc- 
tion in effective cylinder size mentioned earlier, p — 0.45 
corresponds to a cylinder that is approximately half full, 
while for p — 0.3 the layer depth is close to ^D). In 
transient radial segregation is followed by axial, while 
in the other three cases only radial segregation occurs. 
These and other runs suggest that a sufficiently high 
fill level favors permanent radial segregation and at the 
same time suppresses axial segregation. Validation of the 
model, and determining how closely it reproduces exper- 
iment, requires that behavioral trends under a change 
of parameters be reproduced correctly; while experiment 
suggests band suppression at low fill levels [2lj, this ap- 
plies to slurries, and the degree of similarity between dry 
and wet mixtures has yet to be studied systematically. 

The cylinder angular velocity is the parameter most 
readily changed in an experiment. The Q dependence 
(using the values H. — 0.2 and 0.1) for three pairs of runs 
that are otherwise the same has been examined; the run 
pairs are and # J that eventually produce 6 and 3 
bands, #7^ and #i? that produce 12 and 11 bands and, 
lastly, #L and #Cthat both produce 7 bands (a tentative 
result in the case of #(7). A feature common to each of 
these pairs is that Sa is slightly larger at reduced fl. 



The influence of other parameters will be mentioned 
only briefly. Runs #Mand #A^show the effect of remov- 
ing the tangential restoring force, kg = in Eq. ([5]), in 
a sense making the particles less "grain-like" . The runs 
are otherwise similar to and, apart from the func- 
tion used for fy in Eq. now resemble [12]; so do the 
results, in which axial segregation is now accompanied 
by inverse radial segregation. This outcome confirms the 
importance of fg for obtaining the correct radial behav- 
ior. Early in these runs weak normal radial segregation 
occurs, but it becomes inverted when the axial bands 
appear. Different behavior is encountered at a higher fill 
level; run #0, which is related to run #/in a similar way, 
shows normal radial but no axial segregation. Outcomes 
such as these hint at the complex interplay of param- 
eters and call for an extensive coverage, with multiple 
runs, prior to drawing conclusions. 

Another key parameter is b, the relative size of the big 
particles. In run #P, b is reduced to 1.5, and 5 axial 
bands appear, but no measurable radial segregation; b is 
further reduced to 1.3 (the preferred value in [12| ) in run 
=ffQ, resulting in 6 bands. In the opposite direction, it 
is questionable how large b can become without having 
to increase cylinder diameter to avoid serious size effects 
(and, to maintain aspect ratio, also its length), but for 
6 = 2.1 (not shown), depending on fill level, both ra- 
dial and axial segregation can be obtained; experiment 
typically involves & w 3 or greater. 

Making the particles slightly less compressible by dou- 
bling fc„, in run leads to initial radial segregation, 
and eventually three axial bands. A more substantial 
increase in A:„ would, as noted earlier, require a smaller 
6t. Further runs (not shown) similar to #^4, but with 
kn and all the other interaction coefficients increased by 
factors of 10 or 100, and with suitably reduced St, pro- 
duce the same segregation effects, confirming that this 
aspect of the model does not significantly influence the 
behavior. A similar conclusion applies to granular chute 
flow [sHi ■ For a final example of parameter dependence 
(requiring further study), run #5 is similar to #i? except 
that = 7^*" = 10; radial segregation occurs, suggest- 
ing that this effect can be driven by particle size differ- 
ence alone, and eventually two weakly developed axial 
bands appear (the outcome of this long run is less clear 
than usual since a final state had yet to develop) . 



D. Particle motion 

The configurational snapshots recorded during the 
runs enable post-run analysis of the motion of individual 
or selected groups of particles over extended intervals. 
An eventual goal for this kind of analysis is developing a 
capability for relating local organization and dynamics at 
the level of individual particles to the behavior at a scale 
where the collective nature of segregation is exhibited. 
Two examples will be considered here. 

The first is the merging of two bands of big particles. 



FIG. 10: (Color online) Band merging in run #B where small 
particles from the middle band disperse allowing the big par- 
ticle bands to merge; the views show the full system and the 
particles from the bands involved, both before and after the 
event. 



accompanied by the dispersal of the small particles from 
the disappearing middle band. A merge event of this 
type occurs in run #5 (see Fig. ^ between revolutions 
700 and 1200. The upper portion of Fig. [TOl shows three 
views prior to the event, namely the full system, and 
the selected bands of big and small particles (the latter 
also revealing small particles that intrude into the big 
particle band, a detail hidden from the outside). The 
lower portion shows three views at the completion of the 
event, namely the full system, and all the previously se- 
lected particles in their new positions, principally within 
the bands directly involved. Even though there is essen- 
tially no residual small particle core at this stage of the 
run, small particles are able to migrate across the big 
particle bands more readily than the converse; selective 
migration of this kind will be even more apparent in the 
next example. An interesting feature revealed by these 
images is that while the small particles are dispersed in 
both directions, more seem to have traveled to the right 
where the original band of big particles was somewhat 
narrower; while such behavior might not be unexpected 
it is gratifying to see it actually occurring. 

The second example also involves run =ffB, but now for 
monitoring the steady state over a relatively long inter- 
val between revolutions 1230 and 2600. The first view in 
Fig. [TT] is of the entire system at the beginning; the next 
two views show the adjacent, slightly overlapped bands 
of big and small particles selected for tracking. The last 
two views show all the selected particles in their final po- 
sitions (since the band structure is stationary, views of 
the entire system during the interval are practically indis- 
tinguishable); although quite a few small particles (but 
constituting only a tiny fraction of the band population) 
have migrated out of the band in both directions, almost 
no big particles have done likewise, demonstrating that 
they are effectively confined for lack of an escape path. 
Departing particles are replaced by others, but these are 
not shown. This form of analysis can also be carried 
out experimentally by adding tracer particles to existing 



FIG. 11: (Color online) Particle migration and confinement 
over an interval of almost 1400 revolutions during run #_B; 
the views show the full system at the start of the interval 
together with the selected bands of big and small particles, 
and the same sets of big and small particles at the end. 
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FIG. 12: (Color online) Final state of #r; top, bottom, big 
and small particle views. 



bands [22l | , revealing that confinement is indeed a strong 
effect. 



E. Presegregated systems 

An alternative perspective is provided by systems that 
begin in an axially segregated statCja problem also stud- 
ied experimentally with MRI [S^, H^j. Various scenar- 
ios exist, of which preserving the original state and ax- 
ial mixing represent two extremes, while rearrangement 
into some other band pattern is also a possibility. Ex- 
amples of all three types will be considered. Very little 
change occurs over the duration of run #T. Fig. 1 121 shows 
its final state, and from the top and bottom views it is 
apparent that the interface has adopted a curved form. 
The views showing big and small particles separately re- 
veal that comparatively small numbers of each type have 
penetrated into the opposite region; interestingly, in this 
example big particles also appear to be participating in 
the migration. 

Mixing is most easily achieved by setting all the friction 
parameters equal (as in #5, where only a weak axial 
effect was seen). Fig. [T51 shows the outcome of this in 
run #[/; the initial axial segregation vanishes promptly 
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FIG. 13: (Color online) Axial and radial space-time plots for 
#U{nR = 2440, L = 180, L/D = 6) . 




mmjlmnj: 



FIG. 14: (Color online) Axial space-time plot for {ur — 
6110, L = 160, L/D = 4). 



and a certain amount of radial segregation appears. 

The final scenario, the appearance of an altered axial 
band pattern, occurs in run # V (similar to #L but with 
fewer particles), shown in Fig. 1141 The most prominent 
feature here is the splitting of the small particle band 
by a new band of big particles. The process involves 
particle migration in both directions; from the recorded 
configurations (not shown) it appears that small particles 
follow an interior path while big particles tend to prefer 
a path close to the curved cylinder wall; the dynamics of 
this process is yet another aspect of the overall problem 
needing further study. 



F. Surface profiles 

The motion of particles along the upper free surface 
plays a role in the segregation process, although the im- 
portance of this contribution and the mechanism involved 
is open to question. The detailed form of the surface is 




FIG. 15: Averaged surface profiles for big and small particles 
(big and small points) for run #i (D = 40). 

therefore of interest, although it is not readily character- 
ized and, for a specific particle mixture, it will depend on 
fill level, rotation rate and cylinder diameter. The profile 
of the free surface typically varies slightly with position 
along the axis once axial segregation has occurred. For 
these reasons, just one instance is included here. Ex- 
perimental examples, ranging from strongly S-shaped to 
almost linear, appear in [38| . 

Fig. [T51 shows the axially averaged surface profile near 
the end of run ^L, where there are seven bands; big 
and small particles are treated individually, and the re- 
sults represent an average over ten configurations spaced 
six revolutions apart. The profiles are reasonably close 
to linear (curved profiles similar to T?] can be obtained 
at lower fill levels); the big particles appear slightly be- 
low the small because it is the particle centers that are 
monitored. There are deviations adjacent to the curved 
cylinder wall, but here the data binning produces some 
distortion. Whether the slight deviations from linearity 
are important for the segregation process remains to be 
seen. Except at the start, there is no obvious change in 
this profile over the entire run. 



IV. CONCLUSION 

Simulations based on a simple particle-based model of 
granular matter have been used to demonstrate the oc- 
currence of both radial and axial segregation in a rotating 
cylinder containing a two-component granular mixture. 
The earlier inability to reproduce radial segregation cor- 
rectly has been resolved; for the type of model consid- 
ered, this is due to the omission of a component of the 
frictional force that helps resist sliding during collisions. 
The different forms of segregation observed are indica- 
tive of a large, complex parameter space, as is the case 
experimentally. 

In order to learn more about the phase diagram and de- 
termine which parameters dominate in different regimes, 
a systematic study, including multiple realizations due 
to reproducibility issues, is required. The results suggest 
many questions of the "what if" type; of necessity, the 
answers will also have to await further work. The model 
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itself is open to critical examination, in particular, how 
the components of the interparticle force influence the be- 
havior (one example has already been given), as well as 
the relevant ranges of the force parameters. More gen- 
erally, in view of the way static friction is represented, 
how reliable is the model altogether (a question equally 
applicable to similar models appearing elsewhere)? 

In theoretical fluid dynamics, stability analysis [S^ is 
used to determine when symmetry breaking is advanta- 
geous. No analogous theory exists for the granular seg- 
regation phenomena described here, and the question of 
what distinguishes the different segregation modes from 
one another, and from the uniformly mixed state, re- 
mains open; is it, for example, the ability to optimize 
energy dissipation, thereby ensuring the most efficient 



flow? The segregation band coarsening slows with in- 
creasing band width, and the states eventually reached 
have varying numbers of bands and degrees of pattern 
regularity; here the question is whether such patterns 
are stable, or are they slow transients en route to even 
further coarsening? Another issue is the relation between 
the radial and axial segregation processes; are they inde- 
pendent phenomena that just happen to coexist (perhaps 
in competition) in most, but not all situations, or does 
the former actually drive the latter? Finally, no attempt 
has been made to probe the underlying mechanisms be- 
yond the exercises in particle tracking; if simulation can 
be shown to reproduce experiment reliably, then, in view 
of the level of detail provided, it may prove an important 
tool for exploring the complexities of segregation. 
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